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Abstract 

The dynamics of recombination in genetics leads to an interesting nonlinear dif- 
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1 Introduction 



The deterministic limit of the stochastic process of recombination in population genetics 
leads to an interesting nonlinear ODE system that has been studied for a long time. The 
first major advances to understand the classical system are due to H. Geiringer ^1] in the 
1940s (though her formulation was slightly different). A full characterization of solutions 
(though not in explicit form) was later obtained by Lyubich, compare ^2] and references 
given there, and jH] for general background material. 

Motivated by this problem in genetics, a relevant subclass was considered in [Ij, where 
an explicit solution to the resulting large system of nonlinear differential equations was 
constructed. Soon after, this was reformulated as a measure valued differential equation 
and solved for a more general class of state spaces [B]. In this context, the main focus 
was not on the functional analytic properties, but on the solution and the interaction with 
other genetic processes, in particular mutation and selection. 

On the other hand, nonlinear semigroups are rarely known explicitly, and if so, this 
usually rests upon the transformability of the system to a linear one. With this observation 
in mind, it seems worthwhile to also investigate the semigroup aspect separately. To do so, 
we will first consider the most elementary semigroup constituents separately. Among them, 
we will identify mutually commuting ones, which are then used to reconstruct the solution 
to the full recombination equation of [H] in a simpler than the previous way. Furthermore, 
generalizations are possible that are less relevant in biological applications, but illustrate 
that the mathematical analysis may still be pushed considerably. 

In view of the origin of the problem, the formulation is oriented towards the application 
in mathematical biology. Consequently, the formulation is slightly more expository than 
necessary for a purely mathematical audience. In a separate section, an abstract refor- 
mulation of the central observation is added, together with a possible generalization to a 
wider class of nonlinear operators. 

2 Mathematical setting 

If A is a locally compact space (by which we always mean to include the Hausdorff prop- 
erty), we denote by JA(X) the Banach space of finite regular Borel measures, equipped 
with the usual variation norm ||.||, i.e., ||u;|| := |c<j|(A) where \u\ denotes the total variation 
measure of uj. M.{X) is a Banach lattice, and we use Ai + (X) (resp. V(X)) to denote the 
closed convex cone of positive measures (resp. the closed simplex of probability measures). 
For positive measures u, one has ||o;|| = u(X). We will mainly be concerned with the 
situation that A is itself a product space, e.g., A = Ai x A2 with Aj locally compact. Let 
us first recall a well-known result, compare [7] and [HI Fact 1] for details. 
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Fact 1 Let v, v' be two regular Borel measures on the locally compact product space X = 
X\ x X 2 , and let v, v' coincide on all "rectangles" Ei x E 2 where Ei and E 2 run through 
the Borel sets of Xi and X 2 , respectively. Then, v = v' , i.e., v{E) = v'{E) for all Borel 
sets E of X. □ 

The state space X that we need will have a product structure, described on the basis 
of sites or nodes. For later convenience, we use N = {0,1, ... ,n} for the set of nodes, i.e., 
we start counting with here. To node i, we attach the locally compact space Xi, and our 
state space is then 

X = X x X x x . . . x X n (1) 

which is still locally compact. Our dynamics will evolve in the Banach space Ai(X) which 
canonically contains the tensor product M. (X ) ® . . . ® M. (X n ) and also its completion, 
the corresponding projective tensor product. In fact, product measures uj = uj <S> ■ ■ ■ <8> u) n 
with uji G A4(Xi) will play an important role below. 

The main reason for using the above set N of nodes is that we will need ordered 
partitions of iV which are uniquely specified by a set of cuts or crossovers. The possible 
cut positions are at the links between nodes, which are denoted by half- integers, i.e., by 
elements of the set L — {|, |, . . . , ^-j^}- We use Latin indices for nodes and Greek indices 
for links, and the implicit rule will be that a = j s the link between nodes i and i + 1. 

With this notation, the ordered partitions of N are in one-to-one correspondence 
with the subsets of L as follows. If A = {a v . . . ,a p } C L, with < a i+ i, let iV A = 
{A Q , A x , . . . , A p } denote the ordered partition 

A = {0, . . . , [aj },A 1 = { \a{\ [a 2 \ } , ... , A p = { \a p ] ,...,n} 

where [a\ ( \a] ) is the largest integer below a (the smallest above a). So, while A C L, one 
has Ai <Z N for < i < p. In particular, we have N = N and N L = {{0}, . . . , {n}}. With 
this definition, it is clear that N B is a refinement of N A if and only if A C B. Consequently, 
the lattice of ordered partitions of N now corresponds to the Boolean algebra of all subsets 
of the finite set L, denoted by B(L), cf. jH Ch. 1.2]. We prefer this notation to that with 
partitions, as it is easier to deal with. If A C B, we will write B — A for B\A, and A for 
the set L — A. 

Let X, Y be two locally compact spaces with attached measure spaces M. (X) and 
JA(Y). If /: X — > Y is a continuous function and u G M.(X), f.uj := ujo f~ l is an element 
of A4(Y), where f~ 1 {y) '■— {x G X | f(x) = y} means the preimage of y G Y in X, with 
obvious extension to f~ l (B), the preimage of a subset B C Y in X. Due to the continuity 
of /, f~ 1 (B) is a Borel set in X if B is a Borel set in Y. 

From now on, let X = X x . . . x X n and let N and L always denote the set of nodes and 
links as introduced above. Let ■n i : X — > Xi be the canonical projection, which is continuous. 
It induces a mapping from JA(X) to Ai(Xi) by uj i— > tt^uj, where {tx^uj)^) = ujiix^^E)), 
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for any Borel set E C X{. By (slight) abuse of notation, we will use the symbol n i 
also for this induced mapping. It is clear that 7Tj is linear and preserves the norm of 
positive measures. In particular, it maps V{X) to V{Xj) and may then be understood as 
marginalization. 

Likewise, for any index set / C N, one defines a projector 717: j\4(X) — > Ai(Xj) with 
Xf := X i6J ^i- With this notation, Xjy = X. We will frequently also use the abbreviation 
n <a for the projector vr {1 . i[aJ} , and n >a for tt {Nj n} . 

This enables us to introduce a class of nonlinear operators, called recombinators from 
now on. For a G L, we first define an elementary recombinator R a : Ai(X) — > Ai(X) by 
R a {0) = and, if u ^ 0, by 

: = A ((*<«■«>) ® ( f >«' W )) ( 2 ) 
II II 

which is a (partial) product measure. Here and in what follows, we tacitly identify (if 
necessary) a product measure with its unique extension to a regular Borel measure on X, 
which is justified by Fact ^ These elementary recombinators satisfy a number of useful 
properties jHj- In particular, ||/? a (o;)|| < ||o;|| for all uj G A4(X), R a is globally Lipschitz on 
M.(X) with Lipschitz constant < 3, and R a {au) = \a\R a (u) for all a G R and uj G M.(X). 
Moreover, R a maps A4+(X) into itself and preserves the norm of positive measures, so 
that V(X) is mapped into itself, too. Finally, when restricted to Ai+(X), the elementary 
recombinators are idempotents and commute, i.e., R 2 a = R a and R a Rp = RpRa for ah 
a, (5 G L. 

As the next step, let A = {a 1 ,...,a p } C L, with a« < cti+i, and let Na be the 
corresponding ordered partition of N, as explained above. Then, we define the (composite) 
recombinators Ra by 

1 P 

R a (uj) := — (g)(7T Ai .u/), (3) 
ll^ll i= 

again with the continuous extension Ra(0) = 0. Note that R = 1 and Rs a \ = R a in this 
notation. As a direct consequence of the definition, or by a simple induction argument 
based on the properties of the elementary recombinators, one obtains the following result, 
compare f>] for details. 

Proposition 1 Let R A be the {composite) recombinator attached to a subset A of L. 
Then, the following assertions hold. 

1. Il-R^u;) || < ||o;|| for all uj G Ai(X), and R A is globally Lipschitz on Ai(X). 

2. R A is positive homogeneous of degree 1, i.e., R A (auj) = \a\R A {uj) for all a G R and 
uj G M(X). 
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3. R A maps M. + (X) into itself and preserves the norm of positive measures. In partic- 
ular, it maps V(X) into itself. 



4. On M. + (X), one has R A = YlaeA^-a' an d the recombinators satisfy the equat 



ion 



RgRr — Rguh j 



for arbitrary G,H C L. □ 

Note that Part 4, when applied to singleton sets, comprises both the idempotency and the 
commutativity of the elementary recombinators. 



3 Recombination dynamics 

Originally motivated by a problem in biology, compare [3], we are interested in the solutions 
of the nonlinear ODE 

u = := Y^e G {R G -l)(u) (4) 

GcL 

on the Banach space Ai(X), with the restriction g G > 0. In general, even though the 
Picard-Lindelof theorem applies (<P is globally Lipschitz) and the time evolution is given 
by the flow of the ODE, it does not seem possible to write down a closed formula for the 
corresponding semigroup in forward time. 

However, if Qq > only for special G, namely the singleton sets {a} with a 6 L, it was 
shown in j^j that the solution in forward time can be given in closed form, if the initial 
condition is a positive measure. This amounts to constructing an explicit formula for the 
nonlinear flow of this ODE on the positive cone Ai + (X). In what follows, we will give an 
independent derivation of this, and an extension to further solvable cases, independent of 
its biological relevance. Let us start with a general property. 

Proposition 2 The abstract Cauchy problem of the ODE O) has a unique solution. Fur- 
thermore, Ai+(X) is positive invariant under the flow, with the norm of positive measures 
preserved. In particular, V(X) is positive invariant. 

PROOF: Existence and uniqueness of the solution follows from the Picard-Lindelof theorem 
on Banach spaces, see |2J Thm. 7.6 and Remark 7.10], because $ is globally Lipschitz as a 
consequence of Part 1 of Proposition ^ 

Let v e Ai + (X), i.e., v(E) > for all Borel sets E C X. Let E be any Borel subset of 
X such that u(E) = 0. Then 

4>(v)(E) = J2cgRg(v)(E) > 

GcL 
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because each R G {v) is a positive measure and all g G > by assumption, so <L> satisfies the 
positive minimum principle, compare j3J. By a continuity argument, see p. 235, Thm. 16.5 
and Remark 16.6 of [2|, we obtain the invariance of the cone M. + {X) under the flow in 
forward time. 

Finally, observe that <P(v)(X) = for all v G Ai + (X), which is a consequence of Part 
3 of Proposition [T] Let uj q G M. + [X) be the initial condition, and ui t the corresponding 
solution. So, u t G M. + [X) for all t > by the previous argument, hence ||a> t || = u t (X). 
This implies 4:\\u t \\ = <P(uj t )(X) = 0, so ||u;J = ||a> || for all t > 0. □ 



3.1 Time evolution with one recombinator 

Proposition [21 is important because it allows us to concentrate on the positive cone, and 
to study the forward flow, which is a nonlinear positive semigroup. Our first step is to 
analyze the explicit solution, in forward time, of the ODE 

u = 0A (R A -l)(u) (5) 

for an arbitrary, but fixed A C L, with g A > 0. 

Lemma 1 Let v G Ai+(X) and a G [0, 1]. Then, we have 

R A {au + (I - a)R A {u)) = R A (v) . 

So, R A acts linearly on this type of convex combination, and we also get 

(R A -i)(ais + (l-a)R A (v)) = a(R A -l){v). 

PROOF: Since R A is positive homogeneous of degree 1, it is sufficient to prove the claim 
for v G V(X), where \\v\\ =1. If A = {a 1} . . . , a p } C L, and N A = {A , A l7 . . . , A p }, with 
A{ C N, is the corresponding partition of N, we obtain 



R A (av + (l - a)R A (v)) 



v 



(<K Ai .(av + (l-a)R A (u))) = (g) (a7r Ai .v + (1 - a) n Ai .R A (u)) 

i=0 

V 

(air A ..v+{l-a)iz A ..v) = (5?) (vr^.i/) = R A {v) . 



i=0 i=0 

This proves the first claim, while a verification of the second is now straightforward. □ 

Proposition 3 Let the initial condition u for equation be a positive measure. Then, 
the solution for t > is u t = (pf(u Q ) where {<pf | t > 0} is a semigroup of nonlinear 
positive operators. They are explicitly given by 

tpf = e ~ e ^l+ (1 -e- Q ^)R A = R A -e~ e ^(R A -l) 
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and, on Ai + (X), satisfy the equation 

Finally, uj t R A (u ) as t — > oo, for all u G JA+(X) } with the deviation from the limit 
decaying exponentially fast in time. 

PROOF: Let u> be a positive measure, and u t = e~^'w + (l — e~ Q A t ^jR A {uj Q ), Then, we 
have 

<o t = g A e-^ t (R A -l)(uj ). 
On the other hand, since e~ 0At G [0, 1], we can apply Lemma ^ to verify that also 

g A (R A -l)(u t ) = A e-^(R A -l)(u o ) 

which establishes that u t is indeed a solution of the Cauchy problem with initial condition 
uj , while uniqueness follows from Proposition |21 

Let {ipf | t > 0} denote the corresponding forward flow, where iff then is the nonlinear 
operator stated in the proposition. It maps Ai + (X) into itself, preserves the norm of 
positive measures, and satisfies the equation 

by the general properties of the flow. Since the formula for the solution is valid for all 
initial conditions u G Ji4 + (X), the formula for the derivative is correct on A4 + (X). 

The statement about the norm convergence of ui t to the product measure R A (u ) follows 
from the form of (p t by standard arguments. □ 

It is somewhat surprising that such a simple formula for the solution emerges. It admits 
a probabilistic interpretation as follows. If g A is the rate of the recombination process with 
R A , the term e~ B A l (resp. (1 — e~ ^*)) is the probability that recombination has not yet 
taken place (resp. has happened at least once) until time t. This can be substantiated by 
viewing the ODE (J3j) as the deterministic limit of an underlying stochastic process, e.g., in 
the spirit of jTUl Sec. 11]. 

Let us also note that the formula for ipf would formally emerge from expanding 
exp (g A t(R A — 1)) as if R A were a linear idempotent, which it isn't. Nevertheless, the 
result is still correct as a consequence of Lemma [U because linearity on those special 
convex combinations is all that is needed to derive the formula properly. 

Still, the right hand side of (0) is genuinely nonlinear. In general, it is well known, 
compare [SJ p. 91], that T^(f) := / o {p A defines a semigroup of linear operators on the 
dual space of A4(X). One can also define the corresponding generator, but this approach 
does not seem to help in understanding the result of Proposition El 
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However, assume u G M. + (X) and define the two signed measures v x (t) = R A (ui t ) 
and v 2 (t) = (Ra ~ ^)( u t)- Then, one has u t = v^t) — v 2 (t), but also v^it) = Ra^q an d 
v 2 (t) = e~ Q A l (R A — 1)(oj ), by an application of Lemma [TJ Consequently, these measures 
satisfy the differential equations 

i , 1 = and z> 2 = —g A v 2 

which shows that the ODE ((HJ), when restricted to the cone M. + (X), is equivalent to 
a system of two linear ODEs together with special initial conditions. Also, v 1 is the 
equilibrium, and ||^ 2 (i)|| — ► exponentially fast as t — > oo. 

3.2 Abstract reformulation and possible extension 

In view of other applications, it might be worth to extract the underlying structure of 
Proposition |3] as follows. 

Theorem 1 Let K be a closed convex subset of a Banach space B, and let 1Z: K — > K 
be a [nonlinear) Lipschitz map which satisfies 

Tl(ax + (1 - a)TZ{x)) = TZ{x) (6) 

for all a G [0,1] and all x G K . Let g > be arbitrary. 
Then, the (nonlinear) Cauchy problem 

x = g(TZ - l)(x) , x(0) = x G K , (7) 

has the unique solution x(t) = e~ 0t x Q + (1 — e~ 0t )lZ(x o ) for t > 0, and the entire forward 
orbit remains in K . 

PROOF: As above, the uniqueness follows from Lipschitz continuity. The formula for 
the solution is again verified by direct differentiation, using (JHJ). Note that this assumed 
property comprises the linearity on special convex combinations and the idempotency of 
1Z (by setting a = 0). Since K is a closed convex set which is mapped into itself by 71, its 
forward invariance under the flow follows by standard arguments, compare □ 

Several generalizations seem possible at first sight, but most of them break down rather 
quickly as soon as genuine nonlinearity of 71 sets in. First, consider the ansatz 

K(ax + (1 - a)K n (x)) = 7Z(x), (8) 

for some n > 1, and again for all a G [0, 1] and all x G K. This includes the re- 
lation 7Z n+1 (x) = 7Z(x) for x G K, a generalization of the idempotency used previ- 
ously. Here, an explicit solution of the Cauchy problem (J7J) would contain all the terms 
x Q ,7Z(x ), . . . ,7Z n (x ), so that (jHJ) cannot suffice as a substitute for linearity. 
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Alternatively, one might consider the ansatz 

Tl(a x + a x K{x) + ... + a n TZ n (x)) = K(x) (9) 

on K, where the coefficients ai > with a +. . -+a n = 1 define general convex combinations. 
With ctj = 1 (hence ai = for all i ^ j), this comprises W +1 = 1Z on K, hence also 1Z 2 = 1Z, 
so that this ansatz cannot yield an extension of (jOJ. 

Any meaningful generalization of © seems to require a partial linearity of 7Z, namely 
on a subset of the convex combinations used in 0. To construct one possibility, fix n > 2, 
set £ := exp(2-7n/n), and define 

Gt\t) := i(e* + £ fc e^ + £ 2fc e^ + ... + £ ( "- 1)fe e«' 1 " 1 *) (10) 

for k G Z. Due to £ n = 1, this definition is modulo n in the lower index k. Observing 

iy^M = f 1 - if n divides M 
n m=0 [0, otherwise 

and expanding the exponential factors in (|10jh one finds (for < < n) 

G p(t) = ^E^( 1+ ^ + ^ 2(W)+ --- + ^ M(w) ) 

£>0 



5 *,o + E 



jmn—k 



~* (ran — k)\ 

For < k < n, these functions have the following elementary properties. 

2. £Gf (*)=G&( t ) (with G^G^); 
3- ESGi B) (t) = e*. 
Moreover, one has 

Lemma 2 Let n > 2 and < k < n. Then, lim^oo e~ l G^ (t) = -. 
PROOF: If 7 G S 1 , one has 7 = cos(0) + % sin(0) for some G M, hence 

e~*e 7i = e (7 ~ 1} * = e ( cos (^)- 1 )* e * sm(^)t_ 

The absolute value is then e ( c °s(<p)-i)t which tends to as t — > 00, unless cos(</>) = 1. Using 
this term by term in (fTUj) proves the claim. □ 
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Let us assume, as above, that TZ(K) C K and consider the abstract Cauchy problem 
(J7J). Let us also assume, for a fixed n > 2, that lZ n+1 = TZ on K. Define 

n-1 

*>« = e-< (i + (g^w - m n +Y^ G n\m k ) (12) 

*!=i 

so that <£> = 1 and <£> t f ~*°°> ^(72- + 7£ 2 + . . . + TZ n ), as a consequence of Lemma 121 

If we now have [TZ, ip t ] = for all t > 0, one can check, by an explicit calculation, that 

<p t (x Q ) solves the Cauchy problem (|7ji. For n = 2, it is the solution of Proposition 01 

Let us make this approach more concrete for n = 3. Here, one has TZ 3 = TZ, and the 

solution <f t (x Q ) would read 

x(t) = e"'a; + e~* sinh(t) K(x ) + e~*(cosh(t) - 1) ^ 2 (x ). 

The coefficients on the right hand side once more admit a probabilistic interpretation. 
The term e~* is the probability that no "hit" (by TZ) has happened until time t, while 
e~* sinh(t) = |(1 + e*)(l — e"*) (resp. e~*(cosh(t) — 1) = |(1 — e - ') 2 ) is the probability for 
an odd number of hits (resp. an even number > 2) until time t. 

Do such operators TZ exist that are nonlinear? One possibility to construct an example 
is the following. For an idempotent R with R{K) C K, find a map o : R{K) — > R{K) 
with a n = 1 that commutes with R on R{K) (which is nontrivial). If one defines TZ = aR, 
one has TZ n = R and TZ n+1 = TZ on K. 

3.3 Dynamics under compatible recombinators 

Let us go back to the more concrete setting of simple recombination and look at the 
situation of compatible semigroups. If A C L, let 

1(A) := {13 | mm(A) < (3 < max(A)} 

denote the complete stretch of links associated with A. 

Lemma 3 Let A C L be fixed and let a G L be given with a $ 1(A). Let {ipf \ t > 0} 
denote the corresponding positive semigroup according to Proposition^ On A4 + (X) , one 
then has 

rfoR a = R Q o V f 

for allt > 0. 

PROOF: In view of the formula for the operators ipf, it is again sufficient to prove the 
claim on V(X). The extension to M. + (X) then follows from the positive homogeneity of 
the recombinators. 
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Let v be an arbitrary probability measure, and assume that a < min(A) (the case 
a > ma,x(A) is completely analogous). If a G [0, 1], one finds 

R a (au+(1 - a)R A (v)) 
= (a 7i <a .u + (1 - a) n <a .R A (v)) <g> (a n >a .u + (1 - a) n >a .R A {v)) 

a(l - ® (7T>a^) + (7T <a .^) ® (7T >a .i2 A (l/))). 

Observe that, due to our assumptions on a, we have vr <Q ,.z/ = ix <a .R A (v). Using this 
relation twice in the last line above (once in each direction), one obtains 

a*Rjy) + (1 - a) 2 R Au{a} (u) + a(l - a) (R a (u) + R a (R A (u))) 
= aR a (u) + (l-a)R Au{a} (u) = (a 1 + (1 - a) R A )(R a (u)) . 

Since v G M. + (X) was arbitrary, this implies the claim. □ 

Corollary 1 Let A, B C L with 1(A) n 1(B) = 0. Then, on M + (X), 

R B ol P? = vt°R B - 

PROOF: Let ip = al + (1 — a)R A for an arbitrary, but fixed a G [0,1]. Write B = 
{Pi, . . . , j3 r } with j3i < ■ ■ ■ < p r . Then, R B = Y[ i= i Rp by Proposition [TJ4, and a repeated 
application of Lemma El gives 

r r 

*w = (n^) ^ = ••• = ^°(n^) = ^° r b 
i=i ' i=i 

which is valid on M. + (X). This proves the claim because ipf is of the form (al + (l — a)R A ), 
with a G [0,1], for all t > 0. □ 

Theorem 2 Lei A,BcL with 1(A) fl I(-B) = 0. Then, the corresponding semigroups 
commute, i.e., 

<P?°<P? = <P?°<P? 

on M + (X), for all t, s > 0. 

In particular, {(pf o ^ | t, s > 0} defines an Abelian two-parameter semigroup of 
nonlinear positive operators on Ai + (X) . 

PROOF: Let v G V(X). With Corollary [TJ one finds 

((al + (1 - a)R A ) o (61 + (1 - b)R B )) (u) 
= abv + a(l- b)R B (v) + (1 - a) (61 + (1 - b)R B )(R A (u)) 
= abv + a(l- b)R B (v) + 6(1 - a)R A (u) + (1 - a)(l - b)R AuB (u) . 
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Since the last expression is symmetric in (a, A) versus (b, B), the operators in the first line 
commute. 

Since the operators tpf and ipf are of the form used in this argument, for all t, s > 0, 
the claim is true on V(X). By positive homogeneity of the recombinators, it extends to 
allof .M+(X). □ 

This allows to formulate our main result, where we write ipf i for the semigroup attached 
to a set Li C L of links. 

Theorem 3 Let A := Ui<i<r ^ be a su bset of L, with I (LA H I{Lj) = for all i ^ j . 
Then, the semigroups {y?^ | t > 0} mutually commute, and the Cauchy problem 

r 
i=l 

with all q l . > and initial condition u > 0, has the unique solution 

r 

u t = (n^)w- 

Asymptotically, u> t — ^ i2 A (u; ) as t — > oo, t/ie convergence, once again, being exponentially 
fast. 

PROOF: Due to commutativity of the participating semigroups by Theorem |2l one can 
apply Proposition El repeatedly to find 

^ = ( yi {j t rf i ) ° n ^ ) m 

i=l j^j 

= (E^(^-i)of[^)w 

i=l i=l 

r 

= E^(^-i)(^)- 

i=l 

The convergence result towards the product measure R a (uj ) is a multiple application of 
Proposition El together with Part 4 of Proposition ^ □ 

Corollary 2 Under the assumptions of Theorem^ the set 

r 

{[[<pZ\ti>Q} 

i=l 
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forms an Abelian r -parameter semigroup of nonlinear positive operators on M. + (X) . Each 
factor is of the form 

ff- = exp(-^ij 1 + (1 - exp(-^ L i i )) J R Li 
where g L . > is the intensity of the underlying process. □ 

Another obvious consequence is that the forward flow of Theorem which is a one- 
parameter semigroup, is embedded into the r-parameter semigroup of Corollary [21 

3.4 Application to single crossovers and outlook 

In j^j, the biologically most relevant situation was investigated where L was written as the 
disjoint union of all its elements. This resulted in the ODE 

By a different approach, it was shown that the corresponding Cauchy problem with positive 
initial condition u Q has the solution 

<*>t = ^2a G {t)R G (u ) (14) 

GcL 

where the coefficient functions are given by 

a G(t) = n ex p(-^) n i 1 ~ ex p(-^))- 

In our new (and more general) formulation, these coefficients can be seen as the result 
of expanding the product (over a G L) of the commuting semigroups {y?" | t > 0}. It 
is an easy exercise to show that the formula of Proposition |H] then gives an independent 
verification of (j!4)l . 

In this simple combinatorial transformation is possible, namely 

t g ■= £(-i) |H - G| ^, 

hdg 

which admits the inverse R G := Y2hdg^h Mobius inversion, compare ^ Thm. 4.18]. 
It was shown in jB] that the signed measures T G (uj t ) satisfy the linear ODEs 

aeG 
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Since T G (u t ) = b G {t)T G {ui ) with b G (t) = YliH<zG a H^)i u t °f f^l* also admits the repre- 
sentation 

w t = E & gW T g( w o)- 

GcL 

This shows that ()13jl can be transformed to a system of 2' L ' linear ODEs (together with 
a special set of initial conditions), which, a posteriori, provides an explanation for the 
appearance of the "almost linear like" behaviour of the nonlinear flow of the original 
equation (fTHj) . At present, I am not aware of any other example of this kind of "Mobius 
linearization" in the literature. 

A similar observation applies to the situation of Theorem EJ which corroborates the 
intuition that explicit expressions for nonlinear semigroups ought to be related to some 
linear structure. As the above examples show, there are perhaps more possibilities for such 
a connection to be discovered. 

Our above analysis revolved around ordered partitions, which make the treatment 
rather simple and transparent. It is clear, however, that there is no principal reason 
to restrict oneself to this case, and with some extra effort, similar results might also be 
possible for more general partitions. One difficulty here is to find a good formulation for 
the cases where the semigroups (and not just the recombinators) commute. 
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